Finite temperature topological order in 2D topological color codes 
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In this work the topological order at finite temperature in two-dimensional color code is studied. 
The topological entropy is used to measure the behavior of the topological order. Topological order 
in color code arises from the colored string-net structures. By imposing the hard constrained limit 
the exact solution of the entanglement entropy becomes possible. For finite size systems, by rising 
the temperature one type of string-net structures is thermalized and the associative topological 
entropy vanishes. In the thermodynamic limit the underlying topological order is fragile even at 
very low temperatures. Taking first thermodynamic limit and then zero-temperature limit and vice 
versa don't commute, and their difference is related only to the topology of regions. The contribution 
of the colors and symmetry of the model in the topological entropy is also discussed. It is shown 
how the gauge symmetry of the color code underlies the topological entropy. 
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I. INTRODUCTION 

Exotic sates of matter are those that defy the usual 
description in terms of well-known Landau-Ginzburg- 
Wilson paradigm 1 , where an appropriate local order pa- 
rameter characterizes different behaviors of two phases 
on either side of the critical point. These are new phases 
of matter carrying a kind of quantum order called topo- 
logical order— and the transition between various phases 
does not depend on the symmetry breaking mechanism. 
From the experimental point of view, the fractional quan- 
tum Hall liquids^i 5 - exhibits topological order. Differ- 
ent phases of electron liquids carry the same symme- 
try and the phase transition between them is possible. 
Therefore, the Landau theory of classical phase transi- 
tions fails in order to describe these phases. The most 
remarkable properties of such new states are the depen- 
dency of the ground state degeneracy on genus or handles 
of the space, gapped excitations and non-trivial braid- 
ing of excitations. This latter property introduces new 
class of emerging particles with the statistical properties 
that are neither fermions nor bosons. In fact, they are 
anyons. They interact topologically independent of their 
distances much like that of Aharonov-Bohm interaction. 
If by winding of an anyon around another one the wave- 
function picks up an overall phase, the anyons are called 
abelian anyons. But, if evolution of the wavefunction 
is captured by a unitary matrix, the anyons are called 
non-abelian anyons, i.e. they obey non-abelian braiding 
statistics. 

The questions such as what the objects underlying 
the topological phases are and how different phases 
can be classified are still under debate. However, 
some physical mechanisms such as string/membrane-net 
condensation 6 -^, branyons 8 . that are analogous of the par- 
ticle condensation in the symmetry-breaking phases and 
a description in terms of quantum groups^ have been 
introduced. The string-nets are extended nonlocal ob- 
jects and the ground state is a coherent superposition of 
all possible string's configurations appearing at all length 



scale. This physical picture clarifies the topological order 
based on the microscopic degrees of freedom. Emerging 
particles such as fermions or anyons are collective exci- 
tations of strings^ -. 

The fact that in topological ordered phases the ground 
state subspace has a robust degeneracy and the exci- 
tations above the ground state are gapped give them 
the ability of being rigorous quantum memory-ii in the 
sense of error correction. However, these are not suffi- 
cient conditions for being self-correcting code. Perhaps 
thermal noises spoil the self-stability of the cod o 12 ' 13 . 
The construction of fault-tolerant topological quantum 
computation 14 exploits the emerging properties of topo- 
logical ordered phases. The quantum information can be 
stored on the topologically protected subspace being free 
from the decoherence. The robust manipulation of quan- 
tum information is done by braiding of anyon s 15 ' 16 : 17 ! 18 , 
where the unitary gate operations are carried out by 
braiding of anyons. 

The ground state of topological ordered phases is 
highly entangled state and an indicator for topological 
order which is not based on symmetry has been emerged, 
notably the topological entanglement entropy 19,20 . This 
topological quantity appears as subleading correction to 
the entanglement entropy of a convex region. It is a gen- 
eral feature of all gapped phases that the entanglement 
entropy scales with the boundary of the region, the so 
called area law^. However, the appearance of constant 
subleading term is a new feature related to topological 
order. The topological phase can be described as a phase 
of matter for which the low-energy effective theory is a 
topological field theory (TQFT). Topological entangle- 
ment entropy is related to one of the basic parameter 
of the TQFT, the total quantum dimension of emerging 
quasiparticles in the theory. 

The best known model for studying topological order, 
emerging abelian or non-abelian anyons and examining 
its capabilities for topological quantum computations is 
the famous Kitaev's model 2 ^. In the abelian phase, the 
model becomes the well-known toric code with a stabi- 
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lizer structure 2 ^. The stabilizer structure of the code is 
given by a set of star and plaquette operators. Stabilizer 
operators fix a subspace in the Hilbert space of the model 
in which different states are distinguished by means of 
topological numbers. Topological color code is another 
relevant example of stabilizer codes, with enhanced com- 
putational capabilitie s 24 ' 25 . The stabilizer operators are 
only colored plaquette operators. However, both codes 
are topological because the stabilizer operators are local 
and non-detectable errors have non-trivial supports on 
the manifold. Non-trivial string operators stand for en- 
coded logical operators. Topological order in toric code 
and color code are related to different gauge symmetry. 
Indeed, the topological order in the toric code is related 
to the Z2 gauge symmetry, whereas the topological order 
in the color code is related to the Z2 x Z2 gauge symme- 
try. This latter symmetry arises from the contribution of 
the colors in the construction of the code and is respon- 
sible of the string-net structure of the code. However, 
both models have almost the same error threshold when 
one considers the error syndrome measurements^. 

Entanglement properties of both models at zero- 
temperature have been extensively studied. The entan- 
glement entropy of a region with its complement depends 
on the degrees of freedom living on the boundary of the 
region supplemented with a topological term. But, the 
topological term for the color code2£ is just twice than 
the toric code 2 ^ when both lattices are put on a compact 
surface. The difference arises because the excitations in 
the toric code are ascribed to the four superselection sec- 
tors, whereas the contribution of the colors in the color 
code model makes the excitations richer than those of 
the toric code. In fact in the color code the number of 
superselection sectors is sixteen. 

How much topological order is robust against the ther- 
mal fluctuations? This is a very important question from 
the practical point of view since every practical imple- 
mentation of quantum information processing can not 
ignore the effects of thermal fluctuations. Thermal fluc- 
tuations may create excitations (errors) accumulating in 
the system and destroying the quantum information en- 
coded in the states. Original idea of topological order 
which is a zero temperature concept can be extended to 
thermal mixed states 2 ^, where classical limit arises by 
washing out the off diagonal elements of the density ma- 
trix. With such realization of classical topological order, 
it would be tempted to investigate the topological order 
in terms of temperature. It can be used as a measure of 
resilience of a code against the temperature. We would 
like to point out at finite temperature the entanglement 
entropy is no longer a measure of correlations between 
subsystems. Instead, the mutual information does as it 
scales with the boundary of subsystems. The topolog- 
ical order however manifests itself as subleading terms 
of entanglement entropy as well as mutual information. 
As long as we are interested in the topological order, a 
linear combination of entropies or mutual information 30 
of subsystems gives a measure of topological order, since 



in the combination the dependency on bulk and bound- 
ary degrees of freedom are washed out. In the toric code 
each of its underlying loop structure or gauge structure 
contributes exactly half to the topological entanglement 
entropy. This is expected since both loop structures are 
similar in a fashion that they are defined on the original 
lattice and its dual. For a finite size system topological 
order survives even at non-zero temperature o 30 ' 31 . Fur- 
ther increasing of temperature destroys the loop struc- 
tures of the model implying the fragility of the toric 
codei 2 * 3 ^ 3 -. 

In this work we address the above problem, i.e the fate 
of topological order at finite temperature, in the color 
codes. The loop structures of the color code is different 
from that of toric code since they are related to different 
gauge symmetries. We discuss how gauge symmetry af- 
fects the finite temperature properties of the code. We 
attach to each set of plaquettes with the same color an 
energy scale. In the lattice gauge theory these energy 
scales are translated into the chemical potential for cre- 
ating the respective charges. Following the derivation 
of C. Castelnovo, et ai^I, we first impose the hard con- 
strained limit on the string-net structure in o~ z bases and 
allow for the thermalization of the string-net structure in 
a x bases. Then, in order to identify the contribution of 
colored strings in the topological entanglement entropy, 
we fix other loop structures and examine the residual 
topological order. Also in the high temperature limit 
the description in terms of classical topological order is 
recovered. 

The organization of the paper is as follows. In the 
next section the color codes is briefly reviewed. Then 
in Sec piip the thermal density matrix that is needed 
for subsequent arguments is derived. In Sec (|IV|) the en- 
tanglement entropy is derived from the density matrix. 
Then, limiting behavior of the entanglement is discussed 
in Sec((V|. Topological entanglement entropy and its be- 
havior in terms of temperature is the subject of Sec (|VI[) . 
The case of open boundary conditions is discussed in 
Sec dVIip . Sec (|VIII[) is devoted to conclusions. 



II. 2-COLEX AND COLOR CODE: FIXING 
NOTATIONS 

Let us start by a brief introductory on the color code 
model. Consider a two-dimensional trivalent lattice com- 
posed of plaquettes, vertices and links. Such structure is 
shown in Fig[TJ To keep track of vertices, plaquettes 
and links, we use color as a bookkeeping tool. We will 
use three colors: red, green and blue. Then we color 
the plaquettes of the lattice in which two neighboring 
plaquettes don't share in the same color. In this way 
we can also color links so that a c-link (the letter c 
stand for color throughout the paper unless stated an- 
other meaning) connects two c-plaquettes. We call this 
two-dimensional lattice a 2-colex. The lattice can also be 
embedded in higher spatial dimensions, the so-called D- 
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FIG. 1: (color online) A piece of 2-colex with colored hexagons 
as its plaquettes. Lightest to darkest hexagons: green, red and 
blue, respectively (a) The closed green string corresponds to 
the product of blue and red plaquettes. (b) A generalization 
of the closed colored strings, the string-net. (c) Performing a 
a z rotation creates excitations with topological charges (e, \ c ) 
(d) a x rotation creates other types of excitations with charges 
(c,Xe) (e) <r y rotation creates both types of quasiparticles 
at corresponding plaquettes. The composite particle in the 
dashed ellipse is a boson with charge (b, Xb) (f) The composite 
particle in the dashed ellipse is a fermion with charge (b, \g). 



colex^^. Hence, the thermal analysis presented here can 
be extended to higher dimensional cases as in the toric 
code 3 ^. Without loss of generality we suppose all pla- 
quettes are hexagons. Suppose the 2-colex is embedded 
on a compact surface such as torus, i.e imposing periodic 
boundary conditions. To connect the lattice to a physi- 
cal system, we place a qubit at each vertex . With each 
plaquette we attach two operators aimed at reaching the 
stabilizer code. These operators are simply the product 
of a set Pauli operators acting on the qubits around a 

c-plaquette: X c = ILeap °f an< ^ Zc = ILeap°f> wnere 
dP stands for six qubits around a c-plaquette. With this 
definition for plaquette operators we see that all plaque- 
ttes of either color are treated on equal footing. Although 
the colors play no role in the plaquette operators, they 
determine the structure of the coding space. It is simple 
to see that all plaquette operators commute with each 
other since they share either nothing or at even number 
of qubits. The plaquette operators are generators of a 
stabilizer group; thereby there are vectors in a subspace 
of the whole Hilbert space that are common eigenvectors 
of all plaquette operators with eigenvalue +1. This sub- 
space, the coding space, is spanned by some vectors that 
are distinguished by eigenvalues of a set of non-local op- 
erators in which their supports wind the handles of the 
spac o 24 ' 27 . Although three distinct colors introduce the 
code, there is an interplay between color and homology 



of the space. In fact for each homology class only two col- 
ors are independent related to the symmetry of the code. 
As long as we consider closed spaces such as torus, the 
closed non-contractible loops are enough to form bases 
for the coding space. The coding space is called color 
code. From now on, unless it is stated else, we suppose 
the 2-colex spanned by 37V plaquettes, N of each color. 
Notice that on a compact surface like a torus all plaquette 
operators are not independent. In fact they are subject 
to the constraint riggg = Ilbee ^b = IlreTC where 
£1 = X, Z, and Q, B, 1Z are sets of green, blue and red 
plaquettes, respectively. 

The protected subspace is ground state of a many-body 
Hamiltonian that is minus sum of all plaquette operators 
equipped with some coupling constants as follows: 

geg beB ren 

-\ij2 z 3 -^Y, z »- xr *z2 z r> (!) 

geg beB ren 

where \% z 's are coupling constants and each sum runs 
over all corresponding c-plaquettes. The ground states 
spanning the coding space are all vectors in which 
X c \ip) = Z c \ip) — \tp) for all plaquette operators. Any 
violation of this condition amounts to an excited state or 
alternatively as an error. Therefore, the ground state 
is said to be vortex free in the sense of being closed 
string-net condensate. In fact, the ground state is an 
equal weighted superposition of all string-net configura- 
tions. Such configurations can be visualized if we con- 
sider the product of plaquette operators. For example 
consider the product of two neighboring red and blue 
plaquettes. In the product, the Pauli operators of qubits 
shared between two plaquettes square identity and we 
are left with a green string being boundary of two plaque- 
ttes. This closed green string has been shown in FigJTJa), 
which connects green plaquettes. Such closed strings of 
either size and color commute with Hamiltonian. An- 
other interesting extended object underlying the ground 
state structure is string-net that is a collection of col- 
ored strings. A typical string-net has been depicted in 
FigfUb). The associated string-net operator can be the 
product of either a x or a z acting on the qubits it con- 
tains (filled colored circles on the string- net). This string 
operator commutes with all plaquettes since they share 
either nothing or even number of qubits. The appearance 
of such string-net is crucial for the full implementation 
of the Clifford group2i. 

Excitations appear as end points of open colored 
strings. An open string anticommutes with the pla- 
quettes lying at its ends as they share odd number of 
qubits. To have a simple picture of excitations, let us fo- 
cus on some simple cases. Consider a rotation a z applied 
to a certain qubit as in Fig[ljc). As a z anticommutes 
with X r , Xb and X g plaquette operators adjacent to the 
qubit, it will increase the energy of plaquettes by 2A£, 
2AJ, and 2A|, respectively. In this case the excitation on 
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a c-plaquette is revealed by an c-star. Similarly, if we 
perform a a x rotation on a qubit as in Figfljd), the an- 
ticommutation of a x with Z c of neighboring plaquettes 
leads the energy of red, blue and green plaquettes in- 
crease by 2 A!!, 2\ b z and 2Af , respectively. In this case the 
excitation on a c-plaquette is revealed by a c-circle. If a 
rotation a y is performed, as in Fig.fJT];), each plaquette 
contains both above excitations, i.e. the excitations on 
the c-plaquette increase the energy by 2\% + 2\ c z . All 
quasiparticles are by themselves boson. However, they 
may have mutual semionic statistics. For instance, if 
a red particle of Fig.(QJl) go around a green particle of 
Fig. ([Tk) , the wavefunction picks up a minus sign imply- 
ing they are anyons. There are also composite particles 
of two excitations as shown by white dashed ellipses. If 
two excitations differ in both color and type, they form 
a fermion as in Figfljf). Otherwise, they form a boson 
as in Fig[lje). 

All above emerging particles that are the low energy 
excitations of the model can be classified in terms of un- 
derlying gauge group 36 ' 37 . Before that, let us make a 
convention for colors which will be useful for subsequent 
discussions. We refer to colors by a bar operation c that 
transform colors cyclically as f — g, g = b and b = r. The 
elements of the gauge group Z 2 x Z 2 are {e, r, b, g}. Each 
excitation carries a topological charge. The correspond- 
ing topological charge can be labeled by pair (q, x)i where 
q G Z 2 x Z2 and x an hrp of this group. Therefore, there 
are nine bosons labeled by (c, Xe), (e, Xc) and (c, Xc) and 
six fermions (c, Xs) and (c, Xe)- Taking into account the 
vacuum with trivial charge (e, Xe), color code has sixteen 
topological charges or superselection sectors. Now, we 
can label the excitations in FigQ] regarding their topo- 
logical charges. The excitations in FigQJc) and Figfjjd) 
have (e,Xc) and (c, Xe) charges, respectively, while the 
boson in FigQTe) and fermion in Figfjjf ) have (b, Xb) and 
(6, x§) charges, respectively. 



III. THERMAL DENSITY MATRIX OF COLOR 
CODE 

Zero-temperature entanglement of the color code^ 7 - is 
given by writing the ground state of the Hamiltonian in 
Eq.([T]). Let G denotes the group constructed by a set 
of generators of the spin- flip plaquette operators, i.e X c . 
By starting from a polarized state, the ground state will 
be an equal weighted superposition of all elements of the 
stabilizer group as follows. 

i^} = iGr 1/2 £^ o > 0|v| ' (2) 

hEG 

where h is an element of the stabilizer group, \G\ is the 
cardinality of the group, \ V\ is total number of vertices on 
the lattice and a z \0) = |0). Other ground states can be 
constructed through the action of the non-local spin-flip 
operators with non-trivial supports winding the handles 



of the space. They are distinguished by different topo- 
logical numbers. However, all of them have the same en- 
tanglement properties and we consider only one of them 
as above. 

Topological order in the ground state of the color code 
Hamiltonian arises from the condensed string-net struc- 
tures, which are loops without open ends. The ground 
state 1^) has been written in the a z bases. In each term 
of the superposition a collection of spins related to the 
group element h have been flipped. Spins on the cor- 
responding string- net have negative <r z component. In 
this string-net structure Z c = 1 for all plaquettes. There 
are also same string-net structures if we worked in the 
a x bases preserving the X c = 1. This is because the 
model is symmetric upon the exchange of the a z and 
<j x operators on the plaquettes. Unlike the toric code 3 -^, 
both loop structures in the color code exist on the direct 
lattice since for each plaquette we attach two operators. 
The thermal fluctuations can break the string-nets and 
left them with open ends. That how much thermal effects 
are able to destroy the topological order associated to a 
string-net structure depends on the coupling constants in 
the Hamiltonian as well as size of the system. Definitely, 
we calculate the thermal density matrix and the topologi- 
cal entanglement entropy is used as a measure of topolog- 
ical order. Since the ground state is fixed point of all pla- 
quette operators, its underlying topological order arises 
from coherent superposition of two string- net structures, 
each related to one of the bases. We impose a constraint 
in which the string-net structures related to the a z bases 
are preserved. The so called hard constrained limit 29 
can be recasted into the limit X z — > 00. One may ex- 
pect that in this limit the thermal fluctuations can only 
create excitations with topological charges (e,x c )- We 
can describe all excited states as |i/'[7 z ]) — 7 z |'0)j where 
7 Z is any arbitrary product of a z Pauli operators. The 
density matrix p = exp(— (3H) in the bases spanned by 
all states |"0[7 Z ]) is diagonal. However, things get more 
simplified if we work in the original bases /i|0)® '^L Let 
us for simplicity drop the superscript |V|. In these 
bases the density matrix reads as follows 

9 E heG (0\he-^h\0) • [ ] 

Without loss of generality let A z = X z for all plaquettes 
and G x = E g X g , B x = £ fc X fc , Rx = ^ r X r and Z = 
^ c Z c . Then the Hamiltonian proposed in Eq.([T]) can be 
written as 

H = -\ 9 X G X - \ b x B x - \ r x R x - \ 7 Z. (4) 

Since [Z, h] = for all h, then Zh\Q) = 3Nh\Q) and the 
numerator of Eq. © reads as 

(0\he-P H h~h\0) = e 3 ' 3A = Ar (0|e^ G ^^^+' 3A ^^|0).(5) 

Since each plaquette operator squares identity, i.e. X 2 = 
1, each exponential term can be expressed in terms of 
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plaquette operators as follows. 

e PKG x = -Q [cosh{ p K ) + S inh(/3A*)X g ] . (6) 
geG 

The two remaining exponential terms e^ Aa = Sx and e^^ Rx 
have the same expressions and obtained by replacing 
the color index g by b and r. respectively. Every ele- 
ment of the stabilizer group h has a decomposition in 
terms of plaquette operators that it involves as h = 

U g eh X 9 tlbeh X b Ureh Xr > where by c € /i we imply if 
the c-plaquette operator exists in the decomposition of 
h. Let n c be the number of c-plaquettes that exist in the 
decomposition of h. Hence, the only nonzero terms in 
Eq.([5]) are as follows. 

(0|...|0) = (cosh f3X 9 x ) N (cosh (3X b x ) N (cosh (3X r x ) N x 




— {kc-\-k^)N—k c n c -\-kcric-\-k^ri^ _■_ —k g n g —kbnb — k r n r 



(7) 

where k c = — ln(tanh/3A£) and the sum runs over three 
colors c G {r,b,g}. Notice that the operator h not only 
given by the product of green, blue and red plaquettes 
belonging to h but also by other combinations of pla- 
quettes due the periodic boundary conditions. In fact, 
the operator h has the following expressions in terms of 
plaquettes operators. 

g£h beh reh 

= n x > n ^ 

geh beB\h reK\h 

= n *»n** n x r 

geG\h beh reTZ\h 

= II X a II X bU X - ( 8 ) 

geS\h beB\h reh 

Note that in the above expressions by c € C\h we mean 
a c-plaquette in C (C = 1Z,G,B) that is not in h. Now 
we can see that the ambiguity in Eq.((7]) amounts to the 
above expressions for h. The partition function can be 
evaluated in a similar way as follows: 

Z = e 3 ^ N (cosh f3\ x ) N (cosh 0X b x ) N X 

(cosh/3A£)* ^ + 2 e~ {ks+ks)N ^J ■ (9) 

Thus, the density matrix can be recasted into the follow- 
ing form: 

P=W\ £ VT(h)xh\0)(0\hh, (10) 

h,heG 



where 

rfr(h) = 

g— (kc+k^)N — k c n c +kcric+k^n^. _j_ ^ — k g n g —kbnb—k r n r 
J2 C e -{ko+k s )N + i • 

The limiting behavior of the density matrix at zero 
and very high temperatures is compatible with the re- 
spective known results. As temperature tends to zero 
(k c — > 0), the pure density matrix p = J^h h h\0)(0\hh 
is recovered^. At high-temperature limit, the totally 
mixed state p = ^2 h h\0)(0\h yields the classical limit 
of the model upon the hard constrained limit 29 . 

IV. ENTANGLEMENT ENTROPY 

In order to calculate the entanglement entropy we con- 
sider a generic bipartition of the system into subsystems 
A and B. Notice that each bipartition may be com- 
posed of several disconnected regions. Suppose ttia and 
ma stand for the number of respective bipartitions. Let 
Ta and Tb be the number of plaquette operators acting 
solely on A and B, respectively, and let Tab stands for 
the number of plaquette operators acting simultaneously 
on A and B, i.e boundary operators. We focus on the 
entanglement entropy between two partitions A and B 
of the system. To this end, first the reduced density op- 
erator of the one subsystem is evaluated and then the 
entanglement entropy is measured using von Neumann 
entropy. The reduced density matrix of a region, say A, 
is obtained by tracing out p with respect to degrees of 
freedom of the subsystem B. Using the properties of the 
group 2 -, the reduced density matrix reads 

Pa = ]g1 S VT{h)h A \0 A )(A0\hAhA, (11) 

heG,heG A 

where Ga and Gb are subgroups of G which act triv- 
ially on subsystems B and A, respectively. The com- 
plete description of the subgroup Ga is given by a set 
of plaquette operators acting solely on A as well as col- 
lective operators. The latter is a collection of operators 
acting solely on A, but they are not the product of pla- 
quettes in Ga- Suppose A is a simply connected region 
and the set {B\, B2, B mB } presents the disconnected 
components of B. Let ABi be the collection of plaque- 
tte operators acting simultaneously on A and B^. The 
product of c- and c-plaquettes of component Bi with the 
c- and c-plaquettes of the boundary ABi produces a c- 
string acting solely on A. Denoting this string by 7?, its 
expression will be 

7 f= n x*** ( 12 ) 

c,ceBiUABi 

A schematic representation of these collective strings 
is shown in Fig.((2^)- Therefore, for each disconnected 
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region we can realize three collective operators 7^, 7^ 
and 7?. For each disconnected region however only 
two of them are independent due to the local constraint 
lilili = 1 • There is also another constraint on the total 
number of collective operators. It is a global constraint: 
for a given color, say c, the product J\ t 7? can be pro- 
duced by product of all c-and c-plaquettes acting solely 
on A, namely = WcieA-^sX^. Taking into ac- 

count all plaquettes, collective operators and constraints 
on them, the cardinality of the subgroups Ga and Gb 
will be d A = 2 s - A+2mB - 2 and d B = 2 SB+2mA ~ 2 , respec- 
tively. 

The von Neumann entropy as a measure of entan- 
glement between two bipartitions is given by Sa — 
— Tr(p^4 lnp^). We don't use this relation to measure 
the entanglement. Here, we find it useful to instead com- 
pute it using the replica Sa — hmn— >i <9nTr[p^]. Hence, 
the trace of n-th power of the reduced density matrix will 
be^L 




(13) 

In order to bring the above relation to a manageable 
form, we resort to a map between the group elements 
h and the Ising variables^!. Via this map the group 
elements are labeled by Ising variables based on which 
plaquettes and/or collective operators are involved in its 
expression, let 9 C = — 1(+1) be a Ising variables related 
to appearing(not appearing) a c-plaquette in h. Simi- 
larly, 6£ = — 1(+1) related to appearing(not appearing) 
a collective c-string 7? in h. Notice that this map is 
not a one-to-one map. In fact, it is a four-to-one map. 
This latter point arises from the fact that by considering 
three colors, if we reverse the sign of Ising variables re- 
lated to two colors, they represent the same element in 
the group. This point is consistent with Eq.© and will 
be clear in following. The number of c-plaquettes, n c , 
that the element h involved can be expressed in terms of 



Ising variables as 

n c = E^+E(^+^ s J^|^. (14) 

cec i 

Here, Once the color index in n c was fixed, the first 
sum runs over all c-plaquettes in C{Q , B,1Z) and # c 's 
take values of — 1(+1) as defined above accordingly, and 
the E^. is the number of all c-plaquettes in Bi. let 
Yf Bi + t c AB% = E? and since T, C A + £\ E? = N, Eq.fTJ]) 
can be written as 

"-i"^ Sc "^* (15) 

cGC i 

As we stated above, for a disconnected region Bi, the 
strings jT, 7^ and j? are not independent. With such 
dependency, the appearance of 0^0^ in above relation is 
meaningful. To see this, consider the case in which both 
c— and c— strings are present in the group element. This, 
in terms of Ising variables, yields 8f = 6f = —1. But, 
the local constraint on strings implies that the product of 
two strings yield the c— string. Since there is not any c- 
plaquette in the decomposition of a c-string, the number 
of c-plaquettes arising from the collective operators will 
be zero. 

In terms of Ising variables, the Eq. (fT3")) will become 

(7 >. n— 1 n const. 

The coefficient j comes form the 4—1 mapping be- 
tween Ising variables and plaquettes . The sum runs 
over all possible Ising variables subject to a constraint. 
This constraint arises by requiring the non-zero value for 
(0\hi lA h2,A— h n ,A\0) in Eq. (fP3|) . which implies that the 
product of h's must be trivial to give a non-zero value 
for the expectation value. The constraint can be applied 
by the following expression: 



e n 

I r ,Ib,Ig \ i 



s n -m j n e i Q i - i ») 6 [ii ®i e i - t > 



\Z=1 



\l=l 



\l=l 



n 5 ( n e ? - j -)n*(n ^ - * ) n * (n ^ -^y ( w - 1) , (17) 



r \l=l 



b \l=l 



where I c = ±1 subject to I r IbI g = 1 which is imposed 
by the delta function in the second line of above expres- 
sion. Notice that three delta functions in the first line 
are not independent and we could already drop one of 
them. However, we keep them to make the expression 



more symmetric. In fact, once two of them are fixed, the 
third one satisfied. This is just a reinterpretation of the 
local constraint that for a disconnected region only two 
strings arc independent. 

Regarding the map between Ising variables and pla- 
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quettes, we can fully bring the iyr in Eq. ([TO|) into an tjt function reads as follows: 
expression in terms of Ising variables. Using Eq. (|15[) . the 



4Z n 



5{J r J b Jg-l), (18) 



where 

Z = i( e f (fer+fei,+fe s ) + e f (k r -k h -k g ) + 
and J c = ± subject to J r JbJ g = 1 imposed by the 



delta function. Now we are in a position that we can 
rewrite the Eq. (|16p in terms of Ising variables regarding 
to the constraint in Eq. (fTT)) . Also we use the fact that 
nr=iEj^(4) = E {/i}r=i nr=i^({4})- Finally, we 
are left with the following expression. 



\G\ 



n-l 



A 2n Z, 



L E E n E e -^ fl! 



A r =/,. \ / A 6 =/ 6 

n E « 

/ A C ' E ' S = / C 5 s 

„ 1 ^ " e £,(t4s[ef )s e«' )b +^j;s<' )b e<' ) ''e< ,)s +^4s»9f |b 9f )r ) 
V 4 {e< i) - b -'}r = i 



n 



g 2 ^! ^g"g X 



9 {e s }r=i 



(19) 



where 



E 

{4}?=! 

/ 

E 



E E E IPvW-i), 



EEE^/^-i), 

h i 9 

IK, Af = nepep ; c = r> 6,ff. 



The upper limits of sums are just implying the constraint 
in Eq. (|17p . All expressions in the above parentheses can 
be restated as a partition function of a classical Ising 
model. Before that, notice that the index r in the typ- 
ical expression w2^giyl T e~f ^ J r 8 r ^j is mute since the 

sum over configuration of 9 l r 's is done first. Other in- 
dices b, g as well as i are mute. Since the 9 l is an Ising 



variable, it is possible to write it as 



oZ _ _Z_Z+1 



, which 



t s are classical spins. So mapping to the classical Ising 
model provides a tool to have an analytical expression 
for entropy. The constraints A r = +1 and A r = — 1 are 
also satisfied by considering the periodic or antiperiodic 



boundary conditions, respectively. Thus, 

A r =±l 



r m? =1 



1 



2>n (kr,{J r }) 



(20) 



where Z p l a is the partition function of the classical 
Ising model with periodic or antiperiodic boundary con- 
ditions. The partition function of the Ising model can 
be calculated using the transfer matrix 1 as Z = 

Tr[cr^ a 0™=i ^~z]> wnere a x is the usual x-component of 
the Pauli matrices and p/a = 0/1 depending on whether 
we are considering the periodic or antiperiodic boundary 
conditions. The corresponding partition function is 



Z,p/ a = 2™ 



sinh"(^)± (J] J^J cosh "(' 



(21) 



The partition function related to the green and blue pla- 
quettes have similar forms. One needs only replace the 
index r by g or b. The expression in the last parenthesis in 
Eq. (fr9"l) , which is related to the collective strings, can also 
be mapped into a partition function. But, the situation 
is rather tricky. The point is that the terms ©( i )so(') b , 
Q(i)rQ(i)g an d Q(i)bQ(i)r are not all independent. So, it 

is not possible to split the exponential function into three 
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independent terms and then map each term to a partition 
function. However, we can use the idea of mapping and 
transfer matrix. In fact, in this case we use two sets of 
Ising spins. Let choose t, s and it as Ising spins. Setting 
0(090(0* = T l T l+1 , 0(O>"0(O9 = s y+i an d 0(0*0(0'" = 



tt 1 -k l+1 subject to the constraint t l t 1+1 s L s l+1 n l Tr l+1 = 1, 
we hope to calculate the following partition function. No- 
tice that there are a 4-1 mapping between Ising spins and 
colored strings. 



const. , , , 

* * V / k r Tl^rl l + l , k b T l^(l)b I l+l . kg T l ^g ll+l\ 

Z(k r Zr,k b Xlk g Z g i] {4},{J l b },{J l g })= Yl e^ J ^ TT ss + ^ J ^* ). (22) 



First, consider the case with constraints A^ = 1 for c = follows: 
r, 6, g. The transfer matrix will be a 4 x 4 matrix as 



/ e^ b+J 9 g + J - r e - J '» h - jl ! > s+J '' r e jl > h - jl ^- jl *- r e -4 b +4s- J ' r \ 



Ti = 



e -j'b+j'g-j'r e J l b b+j' gg +J l r r ( 



J l b h-J' S -4r 



-J'b+J'g-J'r p J b h+J'g+4r 



-J b b-J'g+4r 



g-j'b-j'g+j'r e J' b b-J sS ~j' r r 



-j'b+j'g-j'r J<b+J'g+J|r 



(23) 



where b = ^E*, g = ^Ef and r = ^Et. The eigen- 
values of the transfer matrix can be easily calculated. Let 
us denote them by & = 4J<£ li , & = 4J<£ 2i , & = 4 46* 
and 4^4i, where 



a* 


= -e b cosh(g - 


f r) 


1 


" b cosh(g - 


-r) 






= -e g cosh(b - 


f r) 


i 

e 

2 


~ s cosh(b - 


t) 






= ^e r cosh(b - 


-g) 


1 

~r 


~ r cosh(b — 


g) 






= -e b cosh(g -+ 


-r)- 


i 

f 2 6 " 


b cosh(g — 


r). 


(24) 



Thus, regarding the constraints A£ = 1 for c = r,b, g 
that are indicated by the superscript ppp in the following 
expression, the partition function reads as follows: 



where the simplification J c — Yi?=i Jc nas been used. 
Considering other constraints that amount to apply- 
ing the antiperiodic boundary conditions on Ising spins, 
the sings of J c in the above expression change corre- 
spondingly. For example, if we were to consider the 
= — A^ = — = 1, the sings of coefficients J g and 
J r become minus. Therefore, other partition functions 
are as follows: 



Z r 7 a = 4^ [J^. - J ff ^ - J r & + . ^ (26) 
Thus, the expression in Eq. (|19p is entirely given in terms 
of partition functions studied here as follows: 
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dn\ 1 



\G\J 4 2 "ZJ 2 S A+2m B 



E E E 

{^>r=i «>r=i {Jj}r=i 



z=l 

(27) 



+ 



+ 



r 



where S, 



yb 



V, 9 



The sum over J^'s can 
be easily taken. Notice J c = ±1. Due to the delta 
function, this leads to J r JbJ g = 1. So, in summation 
the Z 2 x Z 2 symmetry is automatically held. Therefore, 
the factors J c = Yl7=i m the above expression can be 
safely dropped since the mentioned symmetry get simply 
exchanged the terms between brackets. The sums give 
a multiplicative factor (| x 8)™, where the coefficient \ 
arises because of the delta function. It is convenient to 
introduce the notations x c = cosh(-rf) and y c = sinh(^). 
By inserting Eq.(H]), Eq.([25]) and Eq.([26]) into Eq.^i 
we eventually arrive at the following expression. 



1 / cLacLb 



4Z \Zo\G\ 



F, 



(n) 



F, 



(n) 



F, 



A") 



(28) 



where F's are functions of x c , y c and £'s (see Ap- 
pendix (|X A[l ). In particular, the replica trick gives the 
entanglement entropy as follows 



Sa(T) 



( d A d B 
-^-(dF 1 + dF 2 



ln(Zo) 
dF 3 



dF 4 



(29) 



where 9F's are given in Appendix (|X A[) . This relation is 
at the heart of our subsequent discussions. 



V. LIMITING BEHAVIOR OF 
ENTANGLEMENT ENTROPY AND MUTUAL 
INFORMATION 

Eq. ([29]) gives all we need to explore the dependency of 
entanglement entropy on temperature since apart from 
the first term the remaining ones depend on temper- 
ature. First, suppose the size of the system is finite. 



As the temperature goes to zero, all terms that depend 
on temperature vanish and the entanglement entropy of 
ground state is recovered^. The zero-temperature en- 
tropy is S A (T 0) = -ln(^jS-). High-temperature 
limit of the model corresponds to a classical model cap- 
tured through the " hard constrained limit" . As the tem- 
perature tends to infinity, the couplings fc c oo leading 
to x c ~ y c ~ ie fe c/ 2 and & ~ ~ 6h ~ ^ ~ ±e b+r +s. 
At this limit the entanglement entropy acquires a contri- 
bution from the bulk degrees of freedom of region A as 
follows 



S A (T -> oo) = (Yiab + XU - 2toa) In 2. 



(30) 



This is expected in the sense that at high-temperature 
limit the thermal entropy that scales with the volume of 
the region must be reached. This also verifies that the 
extension of von Neumann entropy to finite temperatures 
makes sense, since it contains a bulk contribution (scal- 
ing with the volume of subsystem A) that corresponds 
to the ordinary classical entropy. However, the above 
entropy, despite being at high temperature limit, car- 
ries a constant term. This term 2m a depends only on 
the topology of the region A. This exhibits even at high 
temperature limit the underlying system may carry topo- 
logical order. Here,the classical system is constructed by 
thermalization of a pure density matrix via the high con- 
strained limit 29 . In this way precisely half of the original 
topological order is preserved at the classical limit. Wc 
will refer to this point in next section where the topologi- 
cal entanglement entropy as a measure of the topological 
order is calculated. 

Now, we turn on to take first the thermodynamic limit. 
In this limit, the entanglement entropy behaves as 

S A (T) = - In (^f) + ^(k r + k b + k g ) - ln(4) 



\G\ ) 

(k r + k b +kg 



{dF 1 + dF 2 + dF 3 + dF 4 



(31) 
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(a) 




(b) 



FIG. 2: (color online) (a) A schematic representation of sim- 
ple connected region A and two disjoint regions B\ and Bi. 
For each disjoint region we can realize three colored collective 
strings: blue(dark gray), red(gray) and green(light gray). All 
of them are not independent as in text, (b) A set of biparti- 
tions of the lattice in different manner in which all bulk and 
boundary contribution are dropped through the definition of 
the topological entropy. 



It is simple to check that in the thermodynamic limit 
dF 2 , dF 3 ,dF4 tend to zero, and the above expression be- 
comes 



Sa(T) 



E 



>fcc/2 



(a; c m2; c + y c lny c ) - 



(32) 



To have a more intuitive picture of above expression, let 
us consider the limit of large bipartitions which amount 
to S^,Ep. ^> 1. Finally, we are left with the following 
expression for the entanglement entropy 



Sa(T) 



N- 



(T, AB - 2mA)ln2 



E 



(x c lnx c + y c ln?/ c ) - 



(33) 



If the high-temperature limit is taken, the result coincides 
with Eq. ([3l7)) . However, the result is not consistent with 



the zero-temperature limit. That which limit is taken 
first would affect the entanglement entropy as 



Urn Sa 
lim Sa 



{Yjab - 2TO A )ln2, 

(£ AB -2m A -2m B + 2)ln2. 



(34) 



Both quantities scale with the size of the boundary, i.e. 
the area law holds no matter which limit is taken first. 
However, the subleading terms are different. The differ- 
ence between two quantities is 2(mg — 1), which implies 
that the difference depends only on the topological prop- 
erties of the regions. Like the toric code model^L, the 
topological contribution to the entanglement entropy in 
the color code can be extracted by a single bipartition of 
the model, provided the region B be a multicomponcnt 
region, i.e. raj > 1. 

At zero-temperature the entanglement entropy is sym- 
metric upon the exchange of two subsystems A and B, 
namely Sa{0) = Sb{0)- But, this is no longer true at 
finite temperature. This is because the entanglement en- 
tropy acquires an extensive contribution from the bulk 
degrees of freedom of the region. At finite temperature 
the entanglement entropy is not a measure of quantum 
correlations between subsystems. The relevant quantity 
that drops the bulk dependency and is symmetric be- 
tween two subsystems as well is the so called "mutual 
information" . The mutual information grasps the total 
correlations (quantum and classical) between two subsys- 
tems. In particular, it is a linear combination of entan- 
glement entropy as follows 

I AB (T) = \ (S A (T) + S B (T) - S AUB (T)) . (35) 

Again, the von Numann entropies for subsystems A and 
B as well as whole system AljB are captured by Eq. (|2ll)) . 
Notice that T, A ub = 3JV. The derivation of the mutual 
information is straitforward. However, its full expression 
will be lengthy. Instead of given the full expression, we 
only give the limiting behavior as for the entanglement 
entropy. For finite size systems, the zero-temperature be- 
havior coincides with the entanglement entropy since at 
zero-temperature the system will be in a pure state and 
the entropy Saub{T) vanishes. However, at the high- 
temperature limit, where the classical description holds, 
the mutual information is symmetric between two subsys- 
tems and scales with boundary supplemented with topo- 
logical terms as follows. 

Tab = - (Sab - 2m A - 2m B + 2) ln 2. (36) 

This quantity reveals that in the classical model precisely 
half of the mutual information at zero-temperature sur- 
vives. Apart from the boundary term, the remaining 
terms are topological depending on the topology of the 
subsystems and underlying topological order manifests 
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itself through them. Once again, classical system sup- 
ports half of the topological order of the quantum system 
at zero-temperature. Taking first thermodynamic limit 
and then zero-temperature limit don't commute with the 
opposite case. While both limits scale with the bound- 
ary, the contribution of topological terms arc different. 
The topological contribution to the mutual information 
can be filtered out by taking the difference between two 
following limits: 



lim Iab 
lim 1 A B 

JV->oo,T^0 



(Eab - rriA - ms + 1) In 2, 
(IUb - 2m A -2m B + 2) In 2. 



(37) 



The difference 



AB 



lim 

*o,N- 



Iai 



N- 



lim Tab 

->oo,T^0 



depends only on the topology of the regions as AXab — 
(m,A + rriB — 1) In 2. The advantage of this latter relation, 
in contrast to Eq. (|34p . is that even a simply connected 
bipartition of both subsystems give rises to a nonvan- 
ishing value. So, the mutual information paves the way 
in which we can define a purely topological order quan- 
tity using a simple connected region, in contrast to linear 
combination over different bipartition a 19 ! 20 . 



VI. TOPOLOGICAL ENTANGLEMENT 
ENTROPY 

A. finite values of A % 

Topological entanglement entrop y 19 ! 20 appears as a 
subleading term of the entanglement entropy of a region 
with its complement in the topological phases. It is pro- 
posed that this subleading term can be used as a measure 
of topological order—. In this section we use this mea- 
sure to evaluate the behavior of topological order versus 
temperature in the topological color code. Ground state 
of the color code is appeared as a superposition of two 
underlying string-net structures. Both of them can be re- 
alized on the direct 2-colex lattice (unlike the toric code 
where one of the loop structure is appeared on the dual 
lattice). String-net structures are identified by consider- 
ing the a x or <j z bases. The string-net structure related 
to the latter bases is preserved by imposing the hard con- 
strained limit. But, the string-net related to the former 
bases can be evaporated against the temperature. The 
dependency of the topological entanglement entropy on 
temperature clarifies how much robust the topological 
order related to the string-net structure is against the 
thermal fluctuations. In this subsection we assume that 
all three couplings \% are finite. 

Thermal fluctuations may be so strong that are able to 
break the closed strings and leave them with end points. 



The broken strings as we explained in SecjlTJ carry ex- 
citations on their end points. Imposing the hard con- 
strained limit on the couplings \ c z — > oo in the Hamilto- 
nian in Eq.([T]), the appearance of respective excitations 
is restricted. So, part of topological order related to such 
strings is preserved. As other couplings, A£, are finite, 
one may expect that the thermal fluctuations create the 
star excitations (See Fig[T|) at plaquettes. This means 
that only excitations with charge (e, Xc) are created. As 
we will see the accumulation of these excitations in the 
model destroys the topological order partially. 

As Eg. (f2"9")l suggests, the entanglement entropy scales 
with the boundary as well as bulk degrees of freedom. 
To get ride of boundary and bulk contribution, it is con- 
venient to use a set of bipartitions, and then a linear 
combination of their entanglement entropies unveils the 
topological contribution. This set of bipartitions is shown 
in Fig[2jb). The following linear combination of four bi- 
partitions gives a purely topological contribution to the 
entanglement entropy. 



^topo 



lim (Sia + S 2 a + S 3A - S 4 

R,r — >oo 



(38) 



where each term is given by the von Neumann entangle- 
ment entropy of the corresponding bipartition, R and r 
are the linear size of the subsystems as shown in FigJ^b) . 
Large sizes of the regions are taken to make the size of 
the regions much larger than the correlation length in 
the topological phase. Notice that with the above bi- 
partitions we have the following identifications for the 
number of connected and disconnected regions. 

mi a = m 2 A = m 3A = m 2 B = m 3B = m t B = 1, 

mis = rriiB = 2. (39) 

The following relations between number of plaquettes re- 
lated to different bipartitions also hold. 



E' 



S44 — E 



2A 



J 3A, 



1A 



J l,l 



J 2A 



^3A 



J 4A 



£1,2 = N, 
£2.1 = N, 
£3,1 = N, 
- E5 1 = N. 



(40) 



Notice the ambiguity to the definition of low indices in 
E^ 4 , where j stands for one of the bipartitions in Figj2](b) 



and i has the same meaning as in Eq.([T5]). The full ex- 
pression of Stopo is very lengthy (see Appendix (|X B p ) . It 
is rather hard to see the behavior of the topological en- 
tanglement entropy versus temperature from this equa- 
tion. To make this quantity more clear, let us consider 
some limiting cases. First, we consider the finite size sys- 
tems. Zero-temperature and high-temperature limits are 
as follows: 



T 
T 



(fc c - 
00 (k c 



0) 



OO 



Stopo S cc - 
Stopo S c , 



0. 



-21n2, (41) 
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where S cc — 4 In 2 is the topological entanglement en- 
tropy of color code at zero-temperature. In the zero- 
temperature limit the model coincides with the pure 
ground state. However, at high-temperature limit, as 
it is clear from the above relation, precisely half of the 
topological entropy is removed. This result is obtained 
along the hard constrained limit that has already been 
taken. As we discussed in preceding section, one could 
expect such result. In fact, while in the zero-temperature 
limit the system is fully topological order, at the hight- 
temperature limit only the string-net structures related 
to Z c plaquettes survive. Since both string-net structures 
related to X c and Z c plaquettes have the same contribu- 
tion to the topological order of the ground state, destroy- 
ing one of them give rises to reduction of the topological 
entanglement entropy by half. 

The ground state has topological entropy S cc — lnD 2 , 
where D = 4 is the so called quantum dimension of the 
system. By rising the temperature, the populations of 
open 7 Z string-nets would be favorable since the open 
strings anticommute with plaquette operators X c that 
they share at odd qubits. At the very high temperature 
limit, Sto po (T — > oo) — \nD implies that each underlying 
string- net structure contributes InD to the topological 
entropy of the ground state. 

Now, let us comment on whether the thermalization 
process and taking the classical limit affect the topolog- 
ical sectors of the color code model. The ground state 
of the color code is 4 2g -fold degenerate for the systems 
that live on a manifold with genus g. The 4 2g topological 
sectors arc identified by the eigenvalues of the non-local 
operators winding the handles of the manifold. These 
non-local operators are product of a z operators along the 
winding closed string as S^ z — Ilier c a i ' wnere M an( i 
c stand for homology and color of the respective string 
and T^ is the support of qubits winding the handle. By 
a closed c-string we mean a set of links that connect o 
plaqucttes. So, it commutes with all plaquettes. No- 
tice that for each homology class there are three winding 
strings each of one color. However, because of the inter- 
play between color and homology only two of them are 
independent. Besides, any other closed strings belong- 
ing to the same homology class are equivalent up to a 
deformation. 

Within each sector, the ground state is the equal su- 
perposition of all bases obtained by any given state in 
the sector and applying the plaquette operators as in 
Eq.@. For any ground state, the respective totally 
mixed state that corresponds to the density matrix at 
high-temperature limit is obtained by removing all non- 
diagonal elements of pure density matrix. The value 
of non-local string operators will be preserved by tak- 
ing the high-temperature limit. This latter point im- 
plies that topological sectors are not get mixed through 
the thermalization of the code. In fact, for any mixed 
density matrix, the expectation value of the any closed 
non-winding loops that are product of plaquette opera- 
tors will be +1 in the hard constrained system. Clas- 
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FIG. 3: (color online) Variation of topological entanglement 
entropy in terms of coupling K~S. The horizontal axis is in 
the logarithmic scales. 
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FIG. 4: (color online) Variation of the topological entangle- 
ment entropy in terms of temperature. Different curves corre- 
spond to different sizes of the inner bipartition in Fig[2jb(l)). 
In the limit of high temperatures all topological entropy van- 
ishes. 



sically changing the topological sectors by flipping the 
spins along winding loops are exponentially suppressed 
leading to the so-called broken ergodicity- 2 ^, where the 
phase space is divided into topological sectors. 

What happens if we take the thermodynamic limit 
first? As the size of the system goes to infinity, all 
Hji tends to infinity except one Ef^ related to the in- 
ner region of FigEIb(l)) . Thus, the expression in Ap- 
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pendix (|X Bp gets more simplified as follows 
Stopo (T) —S CC (Q) —> Ar ^°° 

T ^1 1 1 1 



S -x(^ln^ 



fc. 



9 vf 



(42) 



To get more intuitive picture of above relation, we plot 
the variation of S t0 p (T) — S cc (0) in terms of KT, and 
temperature in Fig. and Fig. (J5J) , respectively. For 
simplicity we have set K = k c and X = j . For fi- 
nite size partition, the topological entropy drops as long 
as the temperature increases. However, even at the high- 
temperature limit half of topological entropy survives due 
to the hard constrained limit. As it is clear in Fig [3] at 
KT. = jij, a drop occurs at topological entropy. We use 
this to estimate a characteristic temperature above which 
the topological order disappears. The approximate tem- 
perature reads as follows: 



T, 



Ax- 



drop : 



(43) 



where S' = 3S is the total number of plaquettes in the 
inner region of Figj2fb(l)). Thus, we can conclude that 
the dropping temperature depends on both coupling X x 
and size of the partition. When the size of the parti- 
tion becomes large, the dropping temperature tends to 
smaller values. This is transparent from Figj4] where the 
dropping of the topological entropy versus temperature 
for different sizes of partition has been plotted. Since we 
are considering the large sizes of the partition in the very 
definition of the entanglement entropy, it is natural to 
consider the limit of S ^> 1 in Eq. (|4"2"|) that amounts to 
the following relation: 



Stopo 

(T) -Scc(O) 



N- 



-2 In 2, 



(44) 



which explicitly implies that in the thermodynamic limit 
the topological entropy is fragile at any non-zero tem- 
perature. This is also expected from the dropping tem- 
perature if the bipartitions grow proportional to each 
other. By this we mean that each bipartition scales 
by a coefficient proportional to the size of the system. 
For instance, for the bipartitions in FigJ^b(l)), we can 
write T, C A = 7^-ZV and = 7<iV, where < 7 < 1 and 
7A + 7x + 7a = 1- With this identification, now, we can 
rewrite the dropping temperature in terms of size of the 
system as follows: 



T, 



Ax 



drop 



ln(ViV) 



(45) 



which clearly shows that as the thermodynamic limit is 
reached, the dropping temperature tends to zero. In the 
thermodynamic limit only at zero temperature the topo- 
logical order entirely exists. 



B. finite values of \% but \%, \% -^00 

Thus far, we have considered all A£'s are finite im- 
plying that excitations are allowed in all plaquettes by 
rising the temperature. As we derived in preceding sub- 
section, at the high temperature limit precisely half of the 
topological entanglement entropy vanishes. In fact, the 
string-net structures related to the three colored strings 
are evaporated and excitations in all plaquettes are favor- 
able. In the ground state, the string-net structures con- 
tribute 2 In D to the topological entropy. Now, a question 
arises. How does an individual colored string impact on 
the topological entropy? We give an answer to this ques- 
tion in what follows. To do so, we allow for excitations 
to occur only on c-plaquettes while excitations on other 
c- and c-plaquettes are restricted due to the energy cost. 
This can be done by applying a similar hard constrained 
limit on the couplings A£ and \% in the Hamiltonian of 
Eq.([T]). Let c = g, c = r and c — b. Thus, thermal fluctu- 
ations produce only open green strings carrying topolog- 
ical charges (e,x 9 ). The low lying states of the system 
are IV't'To]) — 7gl'0)j where t| is an arbitrary open green 
string. Within these states, red and blue plaquettes are 
vortex free since [Jf r ,7?] = and [X(,,7?] = 0. 

Taking the limits A£,Afj. — * oo(k r ,ki, — > 0), the topo- 
logical entanglement in Eq. (j4"2"|) behaves as 



Stopo {T) —S cc (0) 



, N- 



_ "9 V y 



(cosh^S^JMcosh^Sf^) 




-(sinh^Sf )1 )ln(smh^Ef A 



As we are eventually interested in large sizes of the bipar- 
titions, i.e. Xf 1 3> 1, any non-zero temperature subsides 
the topological entanglement as follows 



Stopo {T) —S cc (0) 



N- 



ln2. 



(47) 



For finite size partitions and in the high-temperature 
limit, as long as A£ b /T 3> 1, the same reduction in the 
topological entropy occurs. It is obvious from the relation 
in Eq. (|47p that with disruption of one type of strings, say 
green, topological order vanishes partially. This is not 
same as in Eq. (|4"4"]l . where all closed colored strings are 
allowed to evaporate. 

In order to understand the dependency of the topo- 
logical entropy to the colored strings, let us consider the 
non-local winding string operators as order parameters. 
They are products of pauli operators along paths winding 
the torus as S^ x — l~[ier c °i- Consider the ground state 
of the system in the af bases. The excited states are then 
IV'bg]) = IgW- If a non-local string and an open green 
string 7^ cross each other odd times and have different 
colors then {S^' x ,jg} — 0. Otherwise they commute. At 



the ground state level, the expectation values of non-local 
strings are nonzero and independent of any deformation 
of the strings, i.e. T c (0) = {ij;\S^ x \f) = +1(-1), de- 
pending on the topological sector we are analyzing it. 
But, at finite temperature the expectation value is re- 
placed by thermal average as T C (T) = J2{r\ (^i1' x )t- 
The excitations carrying by open green strings are de- 
confined and the expectation value of two non-local red 
(blue) strings on opposite sides of excitations are differ- 
ent. For instance consider an open green string carrying 
star excitations on its ends as shown in FigJS^a). The 
winding blue and red strings cross the green string either 
odd times (solid wavy lines) or nothing (dashed wavy 
lines). The solid lines anticommute with green string 
while dashed lines commute. In the expectation value we 
must take average over all possible cases. So at finite tem- 
perature, the emerging excitations destroy the non-local 
order parameter T r,& (T) ~ 0. However, the expecta- 
tion value of a non-local green string remains finite since 
it commutes with defects. These lead to destroying the 
string-net structure of the ground state through the ther- 
malization. The thermal states still contain closed green 
strings since they commute with defects, i.e. [A^, y£] = 0, 
where is an arbitrary product of a x operators living on 
a closed green string. Notice that closed green strings are 
product of red and blue plaquettes X r and for which 
the expectation values with respect to thermal states are 
+1, i.e. they are stabilized by red and blue plaquettes. 
Thus, one may expect that the topological order in a x 
bases is partially preserved. This is just the message of 
the Eq. (|4"T|) that topological entanglement entropy can be 
reexpressed into St opo = In D + \ In D, where In D is due 
to the string-net structure related to the a z bases and 
\ In D is ascribed to one type of closed colored strings 
(here green) in a x bases survived even at finite tempera- 
ture. 

It is tempting to infer that each colored string con- 
tribute ^hi£> to the topological entropy. However, the 
string structure of the topological order in the ground 
state of the color code is subject to the Z 2 x Z 2 gauge 
symmetry. This is a property of the color code that each 
colored string is a Z 2 gauge degrees of freedom, but all 
strings form a Z 2 x Z 2 gauge structure that is rooted 
in the string-net structure of the model. Nature of the 
topological order in the ground state comes from the fact 
that the ground state is invariant not only under defor- 
mation of strings but also the splitting of a c-string into 
c- and c-strings as shown in FigJSfb). The latter point 
about splitting corresponds to the structure of the gauge 
group of the color code as we explained in Sec. flTf . 

The relation between topological entropy and colored 
strings can be further understood if we soften more cou- 
plings. For example let Ag, \ r x be finite while X b x /T ^> 1. 
Now, the creation of defects in green and red plaque- 
ttes becomes favorable. Although the blue plaquettes 
remain immune against the thermal fluctuations, the av- 
erage value of all colored winding strings vanishes, since 
they anticommute with green or red defects. Thus, all 
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(a) (b) 




(c) 

FIG. 5: (color online) (a) An open green (light gray) string 
carrying star excitations on its ends. The winding blue(dark 
gray) and red(gray) wavy lines on either sides of the excitation 
pick up different values, (b) From top to bottom: a colored 
string, here green, can be either deform or split into red and 
blue strings, (c) A manifestation of a 2-colex with bound- 
ary. Here, we consider a triangular code with three borders: 
blue(dark gray), red(gray) and green(light gray). A basis for 
this code is a string-net in which each string has an end point 
at the border with same color. 



topological order vanishes. This can also be seen from 
the Eq.dU) by letting \ b x /T > l(k b ->• 0). In the limit 
of large bipartitions, all topological order vanishes as fol- 
lows 

Stopo (T) -5 cc (0)^ Jv ^ oo -21n2. (48) 

Once again when two colored strings are allowed to 
evaporate, the topological order in the system subsides. 
Putting all things together, we can conclude that the 
topological entanglement entropy related to the string- 
net structure in the a x bases receives contributions from 
the colored strings. In fact, each colored string con- 
tributes | In D to the topological entropy, but the sym- 
metry of the color code model gives rise to the total con- 
tribution 2 x ilnZ?. 

VII. OPEN BOUNDARY CONDITIONS: 
PLANAR COLORE CODES 

Up to this point, we have only considered the system 
embedded in a closed surface such as torus, i.e. with peri- 
odic boundary conditions. However, the physical system 
must have boundary in which we can confine to a certain 
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piece of space. It is simple to obtain an open surface 
from a closed one by removing some plaquettes of the 
2-colex. Such lattices are called planar color code. For 
example, consider a 2-colex embedded on a sphere. In 
this case the lattice encodes zero qubits, i.e. the ground 
state subspace is spanned by a single vector. If we remove 
a single qubit and its three neighboring links and plaque- 
ttes, the obtained 2-colex will encode a single qubit. In 
fact, three missed plaquettes form three colored borders 
for the lattice, and only strings with the same color of 
border will have end points on that border. An example 
of such bordered lattice is shown in FigJS^c). 

The most important class of the planar color code is 
the so called triangular color code 2 ^. The triangular code 
has three colored borders each of one color (See FigHfc)). 
The logical operators are given by a string-net and its 
deformation in which make a two-dimensional algebra, 
since the code encodes only a single qubit. Such string- 
net is crucial in implementation of the Clifford group. 
The stabilizer group of the code is again given by a set 
of plaquettes. However, all plaquettes are independent 
because of boundary, which affect the number of collec- 
tive operators needed for evaluating of entanglement en- 
tropy. Although a 2-colex on closed and open surfaces 
may represent different properties for the encoding and 
implementation of quantum information processing, they 
represent the same entanglement properties as we will see 
what follows. This is due to the fact that both structures 



have same symmetry and topological order comes from 
the string-net structure of the model. 

The group element h can be produced only in one way 
that is the product of some plaquettes. This amounts to 
consider only one of the terms in Eq. (|7|) and Eq. © . This 
leads to the following simple relation for rjT{h). 

m (h) = e -k g n g -k h n h -k r n r _ ^ 

We can follow the method in Sec (|IV|) to map the contri- 
bution of plaquettes and collective strings in the group 
element into Ising variables. However, we should take 
care about the collective strings. Consider two disjoint 
regions B\ and B2 and connected region A as in FigEJa), 
but, on an open surface. Each disjoint region is sur- 
rounded by three colored strings in which only two of 
them are independent. Observe that two colored strings, 
say red, surrounding regions B\ and B2 are independent, 
in contrast to the closed surface. By this we mean that 
if 7 4 r be a collective string around the disjoint region Bi, 
the rj ; 7J is not the product of blue and green plaque- 
ttes of region A. Implication of this point arises in the 
cardinalities of subgroups Ga and Gb- 

Since the open boundary conditions establish only one 
construction for an group element h in terms of plaque- 
tte operators, the preceding arguments leading to the 
expression in Eq. (|27|) give only the first term as follows: 



1 



1 



\G\J Z«2 S A+2m B 



E E E 



{4>P=i {Jj}"=i { J «}?=i 



1=1 



(z%Xk r ,{4})) A (zto\h,{4}j) A (4 P) (M4))) A n z r p (w,^,fc 3 i]f;{^},{4},{4})|, 

(50) 



where Z\ = e ™ (fc s +fcb+*v)_ j n particular, the replica gives 
the entanglement entropy as 

Sa(T) = In + In(Zx) - j-dF,. (51) 

The expression dF\ is given in Appendix (|X A|) . 

For finite size systems, taking the zero-temperature 
limit gives the entanglement entropy of a single bipar- 
tition, rriA — m-B = 1, in the ground state of the tri- 
angular color code, i.e. S A (0) = (Sab — 2) In 2, where 
the upper index t denotes the triangular code. If we 
were to take first the thermodynamic limit and then 
the zero- temperature limit, the limits don't commute 
with each other. The difference between two limits de- 
pends only on the topology of the regions as AS* = 
limT^o.iV^oo Sa — hnijv->.oo,T->-o S A = 2m B ln2, which 



implies that the topological contribution to the entangle- 
ment entropy can be extracted even by a single biparti- 
tion. 

Topological entanglement entropy is given by Eq. (l551) 
and the respective bipartitions in Figj2^b). We should 
take into account the restriction imposed by bound- 
ary. The cardinalities of bipartitions in FigO(b(l)) are 
d A = 2 SlA+2 and d B = 2 SlB+2 . The corresponding 
quantities for bipartitions in Figj2](b(2)) and Figj2](b(3)) 
are d A = 2 s2 ' 3 - 4 and d B = 2 E2 - 3B+2 . Last bipartition 
in FigJ5fb(4)) yields the cardinalities dA = 2 E4A and 
ds = 2 S4B+4 . This leads to the topological entropy of tri- 
angular color code at zero temperature 2 ^ as S^, c = 4 In 2, 
which is similar with that of the periodic boundary con- 
ditions. 

Finite temperature topological entropy behaves as 
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what we derived for closed surface. For finite size sys- 
tems and at high-temperature limit precisely half of 
the topological entanglement is preserved as Sj opo (T — > 
oo) — 5* c = —2 In 2, which is a consequence of destroy- 
ing of string-net structure in the a x or a z bases. If we 
take first the thermodynamic limit, any non-zero temper- 
ature subsides the topological entropy to half of its value 
at zero temperature. 



VIII. CONCLUSIONS 

In this work the topological order in a class of two- 
dimensional topological stabilizer codes, the so called 
color code, at finite temperature has been addressed. 
Both closed and open surfaces, where the lattice embed- 
ded on, were considered. This is because quasiparticle 
excitations are different for each embedding. The stabi- 
lizer structure of the color code comes from the plaquette 
operators which can be of X and Z Pauli operators. The 
plaquettes are labeled by their color and type. The stabi- 
lizer group is adjusted into a many-body Hamiltonian in 
which the coding space is its ground state subspacc. The 
Hamiltonian is supplemented with energy scales such as 
A^'s and A£'s. The ground state is a string-net conden- 
sate. Closed string-nets are collections of colored strings 
in which no end points are left. Topological order in the 
ground state of the model arises from the coherent super- 
position of two string-net structures that can be visual- 
ized by adopting the string-nets in a z and a x bases. This 
topological order can be characterized by using the topo- 
logical entanglement entropy that is 2 InD where D = 4. 

Considering the limit X c z — > oo, the exact solution of 
the model at finite temperature becomes possible. As 
the temperature is increased, for finite size systems, the 
topological entropy reduced to InD. This implies that 
by thermalization of one of the string-net structures the 
respective topological order is destroyed. Both string-net 
structures have equal contribution to topological entropy. 

The temperature at which the topological entropy is 
dropped is a function of both coupling \ x and size of 
system as in Eg. (1431) . This relation can be used to dis- 
tinguish a length scale for separated defects. It can be re- 
casted into Ne~ 2X */ Tdrop ~ 1. This means that topologi- 
cal order is destroyed when density of defects in the model 
becomes of order unity as the temperature is increased. 
Similar behavior observed in the toric code^i. This in- 
terpretation for density of defects allows us to define the 
respective length scale as ( x = e x */ T . For temperatures 
well below the dropping temperature T <C Td rop , the 
characteristic length scale is large implying the density 
of defects Ae~ 2Ax / T is much less than unity. In this 
case the thermal defects are not able to shave out the 
string-net structure of the model and topological order 
is preserved. As the temperature increases, the density 
of defects becomes of order unity that destroy the topo- 
logical order in the system. In the thermodynamic limit 
and at finite temperature the distance between defects 



is much less than the system size. Thus at this limit 
any nonzero temperature can destroy topological order 
as in Eq. (|4~4")) . The robustness of the string- net struc- 
ture in the hard constrained limit can be understood via 
such identification for distance between defects. In the 
limit A z — > oo the distance between respective defects 
Q z = e Az / T is infinity. So, the respective topological or- 
der is immune even in the high temperature limit. This 
is the reason why in the high constrained system half of 
the topological order is preserved. 

Topological order in the color code arises from coherent 
superposition of string-net structures in a x and a z bases. 
Both structures are similar having the same contribution 
in the topological entropy. By imposing extra conditions 
on the couplings A!^, we can examine the contribution 
of colors and underlying symmetry on the topological 
entropy. If we let for defects to occur only in c-plaquettes, 
i.e. X x , X% — > oo, the thermal states still carry topological 
order in the a x bases even at thermodynamic limit as 
in Eg. (147)) . However, if we soften the latter condition 
and let for defects to occur in both c- and c-plaquettes, 
i.e. only let X x — > oo, as Eg. ([4*5|) suggests topological 
order is entirely destroyed. These observations reveal 
that both colors and symmetry determine the topological 
entropy. In each loop structure, either in a x or <j z bases, 
string-nets are collections of closed colored strings. Each 
colored closed string contributes \ In D to the topological 
entropy. Considering the gauge symmetry of the model, 
the total contributions of closed colored strings in the 
topological entropy will be InD. In the ground state 
both string-net structures in a x and a z bases contribute 
equally yielding the topological entropy as 2 InD. 

All above results also hold for the case of open bound- 
ary conditions. For this case we considered the triangular 
color code with borders. Although lattices embedded in 
closed or open surfaces present different properties from 
the quantum information perspective, they have simi- 
lar topological order. This is because they have simi- 
lar gauge symmetry group and consequently same un- 
derlying string structure that is reflected in topological 
entropy. 

The topological order in both toric code and color 
code models in two-dimensional space is fragile against 
thermal fluctuations in the sense that their underlying 
structures are one-dimensional objects, i.e. the strings. 
This fragility limits their capabilities as self-stability 
codes. The deconfinement of open strings carrying 
excitations can be restricted by coupling the defects to 
bosonic fields^. Toric code model in three-dimension 
have a membrane structure that are robust against 
thermal noises^. The color code model can also be 
generalized to D-dimensional space with branyon and 
brane-net structures. It is instructive to generalize the 
approach presented here to color codes based on the 
D-colexe o 8 ' 34 in the sense that they have different under- 
lying objects as physical mechanism for topological order. 
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The case n = 1 and their derivatives are needed for evalu- 
ation of entanglement entropy. So, the above expressions 
get more simplified as follows: 

F W _ e f(k r +k b +k g )^ F (i) = e f (-k r -k b +k g )^ 
- p ^(k r -k b -k g ) _ p §(-k r +k„-k g ) 



X. APPENDIX 

A. F's expressions and their derivatives 

Fj's for j = 1,2,3,4 stand for four terms between 
the two brackets in Eq.(|27p. By using the partition func- 
tions related to plaquettes and strings that are given in 
Eq.([2"T]). Eq.([2"5]) and Eq.([5BJ, the following relations for 

Fj 1 ^ 's are obtained. 



(52) 



I 

Trivially Z Q = ±(.F\ (1) + F 2 (1) + F 3 (1) + F 3 (1) ), and their 
derivatives at n = 1 are 



F[ n) = (as" + y-f- (x» + y^f- (x» + y^f* J](^ + & + ^ + &) 



F 



(n) 



F 



{n) 



i 

« + y ; l ) E - (i? - tff* (x n g - y-f- IJ(-& - ^ + S + &) 
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dFi = (dF^/dn)^ = Fi 



(i) 



E A e 2 r ( Xr In x r + y r lny r ) + H b A e 2 (x b lnx b + y b In y b ) + Y? A e 2 ( Xg \nx g + y g h\y g ) 



dF 2 = (dF^/dn)^i = FP 



T, r A e 2 ' (av lnx r - y r lny r ) + S^e ^ (x b lnx b - y b lny b ) + Y, 9 A e ^ (x g lnx g + y g lny g ) 



-J2e^ + ^-^ {-Zulntu +6iln6i - 6i lnfo+fo ln&i) 



9F 3 = {dF^/dn)^ = F 



(3) 



S^e 2 " (x r lnx r + y r lny r ) + S^e 2 (ztlnxh - y b \ny b ) + S^e 2 (x s lnx 9 - y g lny g ) 



e ~^> + ^+^ Inf M - £2i ln£ 2l + £ 3 i Info + 6m ln^4i) 



5F 4 = (dF^/dn)^ = F 4 (4) 



S^e 2 " (x r \nx r - y r hxy r ) + S^e 26 (xj, lnxft + y b \ny b ) + T? A e 2 (x g \nx g - y g lny g ) 



E 



e 2 



v ! ~ Ef (f h In 6i - 6i In Ca< - 6w In & + £ 4 i In Ui) 



(53) 



B. Topological Entropy 



Applying Eq.(|29p about each bipartition of Fig.([2p), 
and inserting into Eq. (|38|) . the topological entanglement 



entropy then reads as follows: 
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Stopo{T) — S cc (0) — 

_L £ | e ¥(^ + *(^ + *(^, 4 ) (f (i) l n ^) + 41) l n ^) + l n ^) + m ^)) 

+e *("-^-4<*-^>-*<^ 

^-*(*-^+*("-^^^ 

-^je^^-^H^^iHY^-^^Wta^) +^lnCf +d 2) lnd 2) ) 
+e -¥^-sS.i)-^(^-sS. 1 )+^( J v-sg, 1 ) ( _ c (2) lne (2) +e (2) In^f -efln^ 2) +^ 2) ln^f ) 

+e ¥(W-=S >1 )-^(W-=S, 1 )-^(^- E I, 1 ) ( _^) ^^(2) _ e (2) lne (2) +e (2) lne (2) + ^2) lne (2) } 
+ e -¥(A f -^ 1 ) + ^(^V-^, 1 )-^(^. l)(e (2) ^^(2) _ e (2) lnC (2) _ e (2) In^f +e ( 2 )l ne ( 2 ))| 

_ 1 (e^^-^, 1 )+*^-=S.i)+¥^-=S,x)(f(3) ln^ 3) +^ 3) ln^ 3) +^ 3) ln^ 3) +d 3) l^f) 
4Z [ 

+ e ^(N-^l 1 )-^(N-^ 1 )-^(N-^ l){ _^) lQ ^(3) _ e (3) ln ^(3) +e (3) ^(3) + ^ (3) ^(3)) 
+e -^(iV-E S , 1 ) + i(iV-^ 1 )-M (i v_ sll)(e (3) ^^(3) _ € (3) ln ^(3) _ ^(3) ^(3) + ^(3) ^(3)) J 

+ 1 r e ^ ( ^ i) + ^ (Ar _ sli) + i (Ar _ sli)(e ( 4)ln ^4) +^) ln ^(4) +^(4)^^(4) + ^(4) lne (4) ) 

4z ^ 

+e -¥(^-=;, 1 )-*^-=5,i)+TK^-=S.i)(_ € W ln^ 4) + ^ ln^ 4) - ln^ 4) + ^f^f) 
+e ¥(^-=:, 1 )-*^-=i, 1 )-TF(^-=J.i)(_^) ln£< 4 > ln^ 4) +d 4) ln^ 4) +d 4) InCf) 



where S cc (0) stands for the topological entanglement of References 
the color code at zero-temperature, i.e. S cc (0) — 4 In 2. 
Notice that the upper indices of £'s are related to the 
respective bipartitions in Fig[^b). 
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